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Abstract 

We present some combinatorial results on the stochastic abelian sand- 
pile model. These models are characterized by nondeterministic toppling 
rules. The recurrence checking for the deterministic case can be performed 
using the well known burning test which detects presence of forbidden sub- 
configurations (FSC) in strongly polynomial time. In the stochastic case, 
however, even for Manna's model, which is perhaps the simplest non- 
trivial example, no such procedure is known. In this paper, we address 
the decision problem of the existence of any FSC in a general stochastic 
sandpile. We demonstrate a polynomial time algorithm which, given the 
sandpile graph and toppling rules, decides if there exists an FSC. In the 
event of a positive answer, it generates at least one FSC for the given 
sandpile. Repeated application of the algorithm can be used to find many 
distinct FSCs. We also demonstrate a procedure for creating larger FSCs 
from smaller ones and use this to create FSCs for the Manna's model. 
We hope that the structural analysis of stochastic sandpile we perform in 
this paper, will prove useful in the eventual formulation of a deterministic 
procedure to decide recurrence. 



1 Introduction 

The abelian sandpile model (ASM) is a type of discrete diffusion process defined 
on graphs with a specially designated sink vertex. The sandpile model was intro- 
duced by Bak et. al. [2] as model of self organized criticality in various physical 
systems, including the dynamics of sandpile formation. However, this model 
had non-commutative dynamics which were relatively less amenable to analy- 
sis. Later, Dhar [1] introduced the abelian version of sandpile model and since 
then, it has been an extremely fertile area of research in the statistical physics 
community. It is closely associated with the chip firing game investigated by 
Bjorner, Lovasz and Shor [3J and Tardos flOj. The comprehensive survey article 
by Dhar |6j explores its connection with diverse phenomena such as as stress 
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distribution in earthquakes, size distribution in raindrops, path length distribu- 
tions in loop-erased random walks, for instance. See Kleber [8] for a discussion 
on the proximity of this model with numerous disciplines such as probability 
theory, algorithmics, theory of computing, combinatorics, non-linear dynamics, 
fractals, cellular automata, to name a few. 

In the deterministic abelian sandpile model, "sand particles" are added at 
the vertices of a (multi)graph. A site (vertex) is stable as long as the number of 
particles at the site remains less than its degree. Adding more particles would 
render the site unstable and is accompanied by the unstable site's passing a 
particle along each edge to its neighboring sites. This relaxation process is 
referred to as toppling. One of the sites known as the sink cannot topple. 
To ensure that every relaxation process eventually stabilizes, one needs the 
condition that the sink is reachable from every other site. As the system evolves, 
the sandpile goes through a sequence of configurations. Those which can be 
revisited in any toppling sequence are called recurrent, the remaining ones are 
termed transient. A very basic problem here is to check if a given configuration 
is transient or recurrent. Dhar [1] demonstrated a very elegant method, the 
burning test, which solves the decision problem in time quadratic in graph size. 
The utility of the method is augmented even further as it is an essential tool in 
showing the existence of a bijection between a the set of spanning trees on the 
sandpile graph and the set of recurrent configurations [B] . 

See [5] for an introduction to the stochastic version of the abelian sandpile 
model. In this model, every site has a list of toppling rules and at instability, 
one of these rules is chosen arbitrarily and the toppling proceeds in accordance 
with that. We will define this model formally in the succeeding section. This 
model has a decidedly unique behavior because of the inherent non-deterministic 
nature. However, till now no decision procedure is known which can decide if 
a given configuration is transient. A transient configuration is characterized 
by the presence of forbidden sub-configurations (FSC), which are the minimal 
components that are transient by themselves. Dhar [S] discusses some important 
properties of stochastic sandpile and derives bounds on the weights of recurrent 
configurations. A list of FSCs of Manna's model is also reported, which were 
discovered by empirical investigation. It is still not clear if one can determine 
the complete set of these configurations purely analytically. 

In fact, checking for transience is equivalent to checking for presence of FSCs. 
However, in a general stochastic sandpile, it is not clear if an FSC will exist at 
all. We demonstrate the first polynomial time procedure which decides, for a 
given stochastic sandpile, whether there exists any FSC. Apart from solving 
this decision problem, our proof also sheds light on the structural aspects of 
sandpile graphs. We believe this result will contribute in the eventual resolution 
of the recurrence checking problem. We also demonstrate a simple procedure to 
efficiently generate arbitrarily large FSCs for the Manna's model. 
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2 Abelian sandpile models 



Definition 1. A graph G is an ordered pair {V{G), E{G)) wliere V{G) is called 
the set of vertices and E{G) is a set of 2— subsets of V, possibly with repeated 
elements, the set of edges. 

This is referred to as a multi-graph in literature but we will use graph for 
brevity. The degree of a vertex v V is defined as the number of edges in E 
which contain v. Two vertices v and u are called adjacent (or neighboring) if 
{u, v) G E. A path between two vertices u and v is an ordered sequence of edges 
ei, 62, . . . , efe such that u G ei, u € and for all values of i, n ei+i ^ (j). The 
graph G is connected if there exists a path between any pair of vertices. 

To model an Abelian Sandpile Model, we take a connected graph G with 
a special vertex called sink, denoted s E V . Non-sink vertices in G are called 
ordinary vertices and the set will be denoted by Vo = — {s}. 

Definition 2. The configuration of a sandpile G is a map c : V"o — > N. It will 
be represented as a vector. 

The configuration c tells us the number of sand particles that each of the 
ordinary sites currently contain. The empty configuration is the zero vector. 
The capacity of a site is the maximum number of particles that it can hold. 

Definition 3. :An ordinary node v is said to be unstable in a configuration c if 
c{v) > capacity (w). The configuration c is said to be unstable if any site under 
it is unstable, else it is referred to as stable. 

When a site is unstable it is said to topple, i.e. pass on some of its particles 
to its neighbors. When a site v topples once, it losses capacity(v) particles 
and each neighbor of v acquires a particle for every edge, common with v, that 
appears in the toppling rule. The sink node never topples. We start with the 
empty configuration and keep adding particles one by one on sites of our choice 
and topple when necessary. 

The ASM evolves in time through two modes, particle addition at sites and 
relaxation of unstable sites via toppling. The case of many sites becoming 
unstable simultaneously also poses no complication as the order in which they 
are subsequently relaxed does not effect the final stable configuration that is 
obtained at the end of toppling sequence, hence the prefix abelian. Elementary 
proofs of such confluence properties can be found in the pioneering paper on 
ASMs by Dhar |4j. A toppling sequence is an ordered set of configurations 
where every configuration can be obtained from the previous one by toppling 
some unstable site in it. 

Notation: We write ci > C2 if Vf,Ci(i;) > C2 and ci h C2 if there is a toppling 
sequence which takes ci to C2- Lastly we write, ci C2 if 3c3 > Ci such that 
C3 h C2. We say that a configuration C2 is reachable from ci if ci — ^ C2 and 
unreachable otherwise. In words, one can add particles to certain sites in ci so 
that there exists a toppling sequence leading to C2. Note that being reachable 
is a transitive property, i.e. ci — )■ C2 , C2 — >■ C3 ci — > C3. 
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Theorem 2.1. fJ^.J^) Given any configuration c, there exists a unique sta- 
ble configuration a{c) such that c h a{c), independent of the toppling sequence 
chosen. 

A configuration is called recurrent if it is reachable from any configuration. 
As already mentioned, we say that a configuration Ci is reachable from a config- 
uration Cj if by adding some particles to Cj and subsequently relaxing it, we can 
obtain A configuration is transient if it is not recurrent. The set of recurrent 
configurations is therefore, closed under being reachable. 

Definition 4. A configuration c is recurrent iff Vc' we have c' c. 

Since recurrence persists under particle addition, we have the following prop- 
erty. 

Property 2.1. // configuration c\ < c^, then recurrence of c\ implies that of 

Denote the configuration in which every node v has capacity(?;) — 1 par- 
ticles by 

Cmax • Clearly, given any stable configuration c, one can reach Cmax 
simply by adding the required number of particles at each site. Consequently 
Cmax is a recurrent configuration. Evidently, every configuration reachable from 
Cmax is also recurrent. Checking whether a given configuration is recurrent is 
an important combinatorial problem. For the sandpile defined above, Dhar [1] 
presents a simple and computationally efficient recursive procedure for deciding 
recurrence. This algorithm elucidates the importance of the notion of a for- 
bidden sub-configuration. We defer discussing these till next section. Next we 
introduce the notion of a stochastic abelian sandpile. 

A stochastic sandpile (SASM) is the non-deterministic version of the sand- 
pile defined above. Instead of a single toppling rule, a list of toppling rules is 
associated with each vertex. There is a toppling threshold associated with each 
node, it is same for each toppling rule, but may not be equal to the vertex degree. 
When a vertex becomes unstable, a toppling rule is chosen non-deterministically 
and a toppling, consistent with this chosen rule, occurs. 

Formally speaking, an SASM is a triplet S — (G, T, C). Here G is the graph 
of a deterministic sandpile (with a special sink node) as defined above, T is a 
function which associates each node with its set of toppling rules, and C is the 
capacity function which associates with each node its toppling threshold. So 
T{v) is a collection of multi-sets, such that 

1. Each member of T{v) has at most G{v) elements. 

2. Each member of T(w) contains elements from the set of ordinary neighbors 
of u, i.e. Vo nN(w). 

3. Each edge incident on a vertex v appears in some toppling rule t £ T{v). 

When V becomes unstable, one of the members of T(w), say i, is chosen and 
every element in t gets one particle. In case some node appears multiple times 
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in t, it gets as many particles. Note that if toppling t takes place at v, the 
number of particles going out of v is C{v) and those received by its neighbors 
is \t\. By convention, the difference between these two numbers is the number 
of particles that go to sink. 

In analogy with the deterministic configuration c/ is said to be reach- 

able from a configuration Cj if there exists a strategy of adding particles to con- 
figuration Ci at some sites and choosing toppling rules at each toppling, such 
that the resulting configuration is c/. Consider the special configuration Cmax, 
in which every site v has C{v) — 1 particles. Evidently, this configuration is 
reachable from every other configuration, one just needs to add at every node 
the number of particles needed to make its weight just less then its toppling 
threshold. The configuration Cmax is therefore recurrent. Consequently, all 
configurations reachable from this c„iax are recurrent. The ones that are un- 
reachable are transient. This is the definition of recurrence that we will be 
referring to, throughout this paper. 

Definition 5. A configuration c is recurrent if and only if it is reachable from 

3 Forbidden Sub-Configurations and Recurrence 

Consider the following examples of stochastic sandpile. Both of them are based 
on grid graphs. Each site has stability threshold equal to one, i.e. the sites 
become unstable if they get two particles. The only difference between them is 
the toppling rules associated with them. 

- NS-EW {Manna's Model): This is the North South - East West sandpile. 

It is an n x n grid where each site has a threshold of 1 (that is the site 
topples when there are two particles on it). Upon toppling it either gives 
its vertical (north-south) neighbors a particle each or horizontal (east- 
west) neighbors a particle each. The sites on edges and corners loose 
some particles to sink depending on the number of edges common with 
boundary. 

- NE-SW : In analogy with previous example, this is the North East - South 
West sandpile. Threshold of each site is 1 and when a site topples, either 
a particle is passed to the northern and eastern neighbors each, or to 
southern and western ones. 

The first example, the Manna's model, holds prime importance in being 
the simplest sandpile which captures the essential features of general stochastic 
sandpile. 

Note that we will often observe only a local region in the sandpile graph. 
The configuration of that region is a sub- configuration of the sandpile under 

consideration. It so happens that for any given local region, there exists a set of 
sub-configurations which cannot occur in any recurrent configuration. These are 
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Figure 1: NS-EW Stochastic Sandpile 




Figure 2: NE-SW Stochastic Sandpile 



known as forbidden sub-configurations (FSC). We illustrate this with the aid of 
an example. Consider the Manna's model as defined in the example above. We 
claim that the sub-configuration shown in figure (j3|) consisting of four zeroes is 
forbidden. 

The claim that this configuration cannot arise from Cmax by any sequence 
of particle additions and toppling is easily verified. For a simple proof, consider 
the site that toppled last. Depending on which toppling rule was chosen, either 
the horizontal or the vertical neighbor must have received a particle. It cannot 
happen that starting from a particle at every node, one manages to reach a state 
with no particle on any node. More discussion on this property can be found in 

Given any configuration, deciding the presence of any FSC is equivalent 
to checking if the given configuration is recurrent. Note that one cannot take 
existence of forbidden sub-configurations for granted. There exist stochastic 
sandpile in which all configurations are recurrent and consequently the set of 
FSCs is empty. The second example defined above, of the NS-EW sandpile, is 
one such case. It is a simple exercise to demonstrate that, if all sites topple 
using the same rule then no matter which configuration we start from, the 
empty configuration is reachable from c„iax- So one infers that there are no 
FSCs for this sandpile. In general, whether a given stochastic sandpile will have 
a non-empty set of FSCs is a non-trivial question. In the next section we answer 
this question. This result is later used in generating arbitrarily large FSCs of 








Figure 3: A Forbidden Sub-Configuration of the NS-EW sandpile (Manna's 
model) 
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Manna's model. 



4 A guarantee for the existence of FSCs 

Consider a sandpile (G, T, C) with G{V, E) as the underlying graph. The fol- 
lowing algorithm computes the set of those sites on which some forbidden con- 
figuration can exist. 

Procedure 1 (REDUCE(G, T, C)). The input is a sandpile (G, T, C) and the 
output is the reduced sandpile. 

1. Initialize V to V and F to (j). 

2. Select all nodes v in V such that 3 t e T{v) for which tnS = (j). Call the 
set of such nodes D. If D is empty, then terminate. 

3. Reset V = V - D and F = F U D,. 

4. Repeat step 2 till V becomes fixed. 

5. At termination, the set V contains the sites on which a forbidden con- 
figuration can exist and F the set from which all particles can be flushed 
out. Return the sandpile corresponding to sites in V, such that sites in F 
are deleted from all the toppling rules associated with sites of V and the 
threshold function is restricted to sites in V. 

Based on this procedure, we have the following definitions and properties. 
The algorithm's proof of correctness will follow from these. 

Definition 6 (Irreducible Sandpiles). A sandpile S = (G, T, C) is called irreducible 



if REDUCE(5) = S 

Definition 7 (Sub-sandpile). A sandpile S' (= (G',T',G')) is called a sub- 
sandpile of sandpile S (= (G,T, G)) iff 



wc will also denote this sandpile by S jy. 

Property 4.1. A sub-sandpile is a valid sandpile, i.e. particles are not created 
during toppling and sink is reachable from each site. 

A sub-sandpile is a restriction of the original sandpile to some subset of sites. 
A special case is the sandpile obtained by deleting a site v. Here, v is removed 
from the set of sites and also from all the toppling rules it appears in. Denote 
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Figure 4: An example containing a site which does not belong to any minimal 
irreducible sub-sandpile 



this sandpile by S \v-v The idea behind the notion of irreducibility is that 
while checking for recurrence of a certain configuration, certain sites may not 
contribute to the problem complexity at all. Overlooking these can reduce the 
problem size substantially. The sandpiles which are, structurally speaking, the 
simplest to analyze, are defined below. 

Definition 8 (Minimal Irreducible Sandpiles). A sandpile S — {G,T,C) is 
called minimal irreducible if and only if Vu G V, REDUCE(S' \v~v) = 

That is deleting any site from the sandpile makes its analysis trivial. Fol- 
lowing are some useful properties of irreducible sandpiles. 

Property 4.2. For a sandpile S containing a site v, there may not he any 

minimal irreducible sub-sandpile containing v. 

As an example, consider the sandpile in figure Here, all the sites are 
shown with the list of toppling rules associated with them. The sites with single 
headed arrows have unit threshold, and those with double sided arrows have 
threshold 2. It is clear that no matter which site we remove, the REDUCE 
procedure will leave only one or both the square blocks at the end, which means 
there doesn't exist any minimal irreducible sub-sandpile which contains v. 

Property 4.3. The empty configuration on any irreducible sandpile is forbid- 
den. 

Proof of Property 14. 31 : The proof of this is by contradiction, it is a simple 
extension of the idea behind demonstrating the smallest FSC of Manna's model 
in previous section. Consider the last site which toppled. No matter which 
toppling rule was used, there must be a site in the sandpile which should have 
received particle as a consequence (since the sandpile is irreducible). Hence the 
configuration containing all zeroes cannot be achieved via toppling and particle 
additions to Cmax- O 

The structure of the sandpile (in context of reducibility) and its weight (total 
particle count) behave in opposite directions with regards to the recurrence of 
a configuration. Property 14.31 poses a lower limit on the structural complexity 
of a sandpile when the weight is zero, to ensure transience. The next property 
is complementary to the previous one, in that it poses an upper limit on the 
weight for a given structural complexity. 
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Property 4.4. The configuration with one particle on any site of a minimal 
irreducible sandpile is recurrent. 

To prove this, we will show that this configuration is obtainable from a 
configuration in which all sites are unstable and is consequently reachable from 
Cmax- Let S — (G, T, C) be the irreducible sandpile and v be any given vertex. 
Run the algorithm REDUCE on the sub-sandpile S\v-v Since S is minimal 
irreducible, the algorithm will terminate with a null set. Denote the set removed 
from S in the i^*^ iteration by Di. Then Di& have the following property. 

Property 4.5. \/vj e Di, 3t G T{vj) such that tf] Dk = (f) for all k > i. 

This means that each of the vertices removed in z*'* iteration have at least 
one toppling rule which affects only those vertices which have been removed in 
preceding iterations. The validity of this property follows from the definition of 
algorithm trivially. 

Proof of Property 14.41 : For every vertex Vi define the sets Di^i, Di,2, 
. . . , Di fc. , where ki is the number of iterations taken by the algorithm to reduce 
the sub-sandpile to null, as above. Now consider the configuration c in which 
only one site, say Vm, has one particle. All other sites are vacant. 

c{vj) = S^j 

here, c{vj) denotes the number of particles at vj. Si,j is the Kroneckor delta 
function, its value is 1 when i — j and zero in all other cases. We will now show 
that, 

1. c is obtainable from a configuration c' in which all sites in Dm.i can have 
arbitrarily large heights 

2. Any c' in which sites in Dm.i, £'m,2, Dm.k have arbitrarily large 
heights is obtainable from a configuration c" in which all sites in Dm,i, 
Dm,2, • ■ • , Dyn,k, DmM+i havc arbitrarily large heights 

These two properties would imply that c is obtainable from a configuration 
in which all sites have arbitrarily large heights. We first prove the first part, ci 
has one particle at «,„ and zero elsewhere. Without loss of generality assume 
that m = 1. 

ci{vj) = 5i,j 

By minimal irreducibility of S, \ Di i |> 1. Choose any member of Di.i, say 
V2- The definition of Di i implies 3t £ T{v2) such that vi E t . Therefore, using 
this toppling rule, ci is reachable from the following C2, 



9 



On similar lines, let V3 be an element of -D2,ij then C2 is reachable from the 
following C3, 

C3{v,) = {C{v2)-l).S2,j+C{vs).Ss,j 

In general, Cfe is reachable from the following Ck+i, 

Ck+i{vj) = ^c{vi).6ij + {c{vk) - i)-Skj + C{vk+i).Sk+i,j 

i<k 

here Vk+i € Dk,i. The finiteness of G implies that some Vk is revisited at 
some step m. So, 

Cm{'"j) = ^{C{Vi) - l).6i^j + C{Vk).Sk+l,j 

i<'m 

Here, Vm is the same as Vk that is the loop is formed at vertex Vm- From the 
corresponding configuration Cm, repeat the whole procedure again, except that 
now at each step, choose from those members of the Dj.is under consideration, 
which have been chosen the least number of times. In maxi \ Di^i \ iterations, 
every vertex in every Di^i will have at least C{v) — 1 particles. By repeating 
the whole process Z.rriMXi \ Di^i \ times, every node in every Di^i will become 
unstable. We have assumed that for all nodes L{v) > 1, which is a realistic 
assumption, since every site in the sandpile has non-zero capacity. 

L{v) > 2 ^ 3.(L(w) - 1) > L{v) 

This completes the proof of the first part. The proof for second part follows 
the same idea. Lot c be a configuration in which nodes from -Dm.i, -Dm, 2, • • • , 
Dm,i have arbitrarily large heights. Then c is obtainable from the following ci , 

\/Vi e Dm,i+i^t G T{vi) such that 
t n D^^i+j =(/, Vj > 1 

choose one such t for every site Vi in Dm,i+i and denote it U. Let the size 
of Dm,i+i be r. Then the configuration is obtainable from a configuration C in 
which all sites in _D,„,;+i are unstable and heights of sites in Dm,i, Djn,2, ■ ■ ■ , 
Dm,i is reduced by at most r.{maXy^^D,n,i+i^i'"i))- Thus the second part is also 
proved. So, for every minimal irreducible sandpile, any configuration with just 
one particle on any site is recurrent. □ 

Definition 9 (Weight Maximal Transient Configurations). A configuration c 

over a sandpile S is weight maximal transient if it is transient and adding any 
single particle to it at any site makes it recurrent. 
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Corollary 4.1. The empty configuration over a minimal irreducible sandpile is 
weight maximal transient. 

Remark: The correctness of procedure REDUCE can be proved using essen- 
tially the same principles, which imply the properties elucidated above. The 
procedure terminates with a sub-sandpile such that for every site in it, every 
toppling rule effects some other site in the sub-sandpile. The empty configu- 
ration on such a sandpile is transient. This can be shown if one considers the 
last toppling before the empty configuration was attained. No matter which 
site toppled using whichever toppling rule, some other site must have received a 
particle. Hence, empty configuration is not reachable. The procedure therefore 
produces a set of sites, on which the empty configuration is transient. 

5 Generating FSCs for Manna's model 

In this section we will outline a scheme to generate arbitrarily large FSCs for 
the Manna's model (the NS-EW stochastic sandpile defined previously). Dhar 
[5] discusses the FSCs of Manna's model in detail and mentions some of them. 
These FSCs were derived empirically by running extensive simulations. We will 
outline a scheme which efficiently generates most of these. 

Property 5.1. If Hi and H2 are irreducible sub-sandpiles, then 

1. Hi U H2 is irreducible 

2. Any configuration over Hi U H2 with just one particle is forbidden. 

3. ( Gluing Property of forbidden configurations ) If the following conditions 
hold 

(a) ci is a forbidden configuration over Hi and C2 forbidden over H2 

(b) ci and C2 agree over the intersection of Hi and H2 that is Ci{v) = 
C2{v)\/v e V(Hi)r\V{H2) 

denote by ci Vc2 the configuration on V{Hi) UV{H2) which agrees with ci 
on V{Hi) and C2 on V{H2). Then, the configuration obtained by adding 
an extra particle to ci V C2 at any G V{Hi) fl V{H2) such that $t G 
T(vk) with t — {vi} is also forbidden. Basically, we are adding a particle 
at a site in V{Hi) n V{H2) which is not connected to any other site in 
ViHi)nViH2). 

Proof of Property 15.11 : The proves are as follows 

1. Since both Hi and H2 are irreducible, any site in H1UH2 belongs to at least 
one of the component sandpiles and hence has at least one member from 
every toppling rule in the composite sandpile. Therefore, it is irreducible. 

2. Empty configurations over each of the sub-sandpile is forbidden (property 
14. 3[) . The given configuration over the union contains at least one of these 
as a sub-configuration. Hence it is also forbidden. 
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Figure 5: A Forbidden Sub-Configuration of the Manna's model 
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Figure 6: A Forbidden Sub-Configuration of the Manna's model 



3. In all configurations which can generate a configuration satisfying property 
I5.1l bv one toppling, the particle must come from some site in either }i\ or 
_ff2- AH these configurations, therefore must contain either c\ or C2 as a 
sub-configuration (over ll\ or Hi-, as the case may be); both of which are 
forbidden. Hence, the configuration under consideration is also forbidden. 

□ 

Using this property, we will now proceed with outlining our procedure for 
generating FSCs for the Manna's model. Recall that this is just the grid based 
sandpile where each site has a unit threshold and can topple vertically (the 
vertical neighbors get a particle each) or horizontally (the horizontal neighbors 
get a particle) each. The smallest FSC is shown in figure ([3]). We will use 
two copies of this to create a larger FSCs using the property 15.11 See figure 
[SI The two smaller configurations are arranged so that they have a common 
comer site. An additional particle at this site does not make the configuration 
recurrent. It is easy to check that these configurations and there union satisfy 
all the preconditions required by property 15.11 This gives us our first FSC, 
which is not minimal. Repeated application of this rule, gives us a yet larger 
FSC. See figure H 
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6 Future Work 



The main open question is that of coming up with a polynomial time recurrence 
checking algorithm. A randomized procedure with some performance guaran- 
tee will also be a highly significant. However we believe that the logical next 
step from here is coming up with an efficient scheme to generate all the for- 
bidden configurations. There is exactly one FSC in Dhar's list [5] which our 
procedure can not generate. It is not known if there are other configurations 
like this. Improving the scheme we have proposed to account for these missing 
configurations is going to be a crucial step in tackling the original problem. 
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